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■ Abstract. I review existing methods for determining mass-loss rates of red 

O | giants and red supergiants based on infrared data. The simplest method is 

D , based on models for the absorption and emission by dust which forms in the 

f"*^ ■ dense outflows from these cool stars. I discuss the parameters and assumptions 

upon which the method relies, review relationships between the mass-loss rate 
and infrared colours or far-infrared flux density, and propose a new formula for 
the mass-loss rate as a function of the visual extinction. I also briefly discuss 
the use of atomic and molecular transitions at infrared wavelengths. 
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1. Introduction 



q , Early spectroscopic observations of cool giant stars revealed violet-displaced 

absorption in the optical line profiles of strong electronic transitions of singly- 
ionized and n eutral atoms, and mass- loss r ates of M ~ 1CP 8 to 10~ 6 M© yr _1 
were inferred ( Deutscb]|l956l ; lReimerslll975l ). More evidence of such winds came 
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in the form of radio emission lines from abundant molec ules, in particular the am- 



plifie d stimulated emission from hydroxyl at A = 18 cm (jElitzur. Goldreich. &; Scoville 



19761 ) . Thermal line emission from carbon monoxide at mm wavelengths enabled 



measurements of the wind speed and mass-loss rate fo r carbon stars as well as 
stars with oxygen-rich envelopes ( Knapp Sz Morrislll985l ). Methods for determin- 



ing the mass-loss rate using radio techniques are described by Fredrik Schoier 
elsewhere in these proceedings. 

Infrared surveys revealed stars which are heavily reddened by circumstellar 
dust. The association between these dust-enshrouded stars and masers was 
quickly mad e, and the accompany ing dust emission became a very useful tracer 
of mass loss ( Gehrz Sz Wool! 197ll ). It also explained the pumping of the masers 



through infrared emission, but most importantly it provided a mechanism for 
driving the wind. Strong pulsation of the cool photospheres of these stars acts 
as a piston to increase the s caleheight o f the molecular atmosphere, facilitating 
the condensation of grains ( Jural fl986l ) . The continuum opacity of the grains 



allows for efficient transfer of momentum from the stellar radiation field onto 
the dust grains. If the density is high enough then the dust and gas are coupled, 
mostly via grain-H2 collisions. Hi gh mass- loss rates were estimated for these 
stars, M ~ 10" 6 to 10" 4 M yr" 1 (|Qlnon et al1ll984l b 



The Spitzer Space Telescope has made it possible to detect the infrared 
emission from circumstellar dust around red giants and red supergiants (RSGs) 
in many Local Group galaxies for the first time. This enables the study of mass 
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loss and dust production in a rich variety of environment s inaccessible before , 
for example in the Wolf-Lundmark-Melotte dwarf galaxy ( Jackson et al. 20071 ) 



but also galactic g lobular clusters such as the extremely metal-poor Messier 15 
(|Bover et alJl200d i. whilst more thorough and detailed studies are possible in the 
Magellanic Clouds. As technological progress inflates the observable Universe, 
much of the groundbraking work is done in uncharted territory and at the limit 
of sensitivity, running into many of the types of problems that were faced in the 
1970s in studies of nearby stars. It is thus essential to understand the methods 
that are applied to measure mass-loss rates from these infrared data, and in 
particular to be aware of their assumptions and limitations. 



2. Dust-driven wind model as the basis for measuring mass-loss rates 

Circumstellar dust grains absorb stellar light mainly at optical (and ultraviolet) 
wavelengths, and re-emit it mainly at infrared wavelengths. The shape of the 
observed spectral energy distribution (SED) depends on the optical properties of 
the grains and on the optical depth of the envelope. For a spherically symmetric 
geometry the integral under the SED yields the bolometric luminosity, provided 
that the distance is known. The SED does not allow a direct measurement to 
be made of the mass-loss rate. Under the assumption of radiative equilibrium 
between the capture of photons and isotropic emission by the heated grain, 
and applying the continuity equati on, one obtains a cru de relationship that 



summarises the problem quite well (jlvezic Sz Elitzurl ll995): 



r °c t=, (1) 

v exp VL 

where we notice that although the optical depth, r, is proportional to the 
(gas+dust) mass-loss rate, M, it also depends on the dust:gas mass ratio, ip, 
the expansion velocity of the wind, v exp , as well as the luminosity, L. 

Simple radiation-driven dust wind theory predicts how the wind speed 
should depend on the luminosity and dust:gas ratio. This can be used to compute 
the expected wind speed for stars for which we have no direct measurement of 
it. The momentum equation relates the motion of the matter and photon fluids: 

M«exp oc tL, (2) 

where the optical depth properly accounts for the scattering of photons off the 
circumstellar grains. Hence, combination with Eq. (1) yields: 

4, 



HL. (3) 

It seems deceptively reasonable that the dust:gas ratio depends on metallicity. 
Unfortunately, this is a difficult parameter to measure directly. 

First eviden ce for the metallicit y dependence of the wind speed of red giants 
was presented b y Wood et al.1 ( 19921 ) who detected six OH/IR stars in the LMC. 



Marshall et al. I (120011) enlarged this sample, and from a comparison with the 



wind speeds and luminosities of OH/IR stars in the galactic centre they could 
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confirm Eq. (3) and indicate that the dust:gas ratio is linearly proportional to 
metallicity: 

i/) oc Z. (4) 

In the absence of direct measurements of wind speed and dust:gas ratio, it is 
therefore recommended to abide by the following scaling relation: 

— oc Z-°- 5 L - 75 , (5) 
r 

The metallicity dependence of dust-driven winds is reviewed in Ivan Loonl (|200fih . 



3. Model computations of the infrared spectral energy distribution 

The optical depth of the dust envelope can be derived from the observed SED 
by comparison with synthetic SEDs. The latter are produced by computing 
simultaneously the radiation transfer through the dust envelope, and the ther- 
mal balance between the grain heating through irradiation and cooling through 
reradiation. The prin ciples underly i ng these computat i ons w ere fir st outl i ned fo r 



reradiation. ine prin ciples underly i ng tnese computat i ons w ere nr st outl i ned lo r 
the non-grey case by iLeung owan- Robinson! l|l980l ) and lYorkd ([1980). 



Three commonly used codes that solve the equations of radiation transfer and 
thermal equil i brium were developed in dependently in the 1990s: the co de of 
Groenewegenl (ll993h. the DUSTY c ode (jlvezic. Elitzur. Nenkoval Il999 ), and 



the modust code ( Bouwmanll200ll ). Of these, only DUSTY is publicly available 



fortunately it is also well documented and relatively easy to use. The latest 
generation of codes are based on Monte Carlo techni ques for tracing the energ y 



a,cmg tne energ y 

packets on their journey through the dust envelope (jBiorkman Woodll200ll ). 
culminating in t he publicly available photo- ioniza tion /dust radiation transfer 
code MOCASSIN (jErcolano. Barlow, fc Storevll2005l ). 



As an illustration, the dusty code is used to generate a series of model SEDs 
with a visual optical depth between ry = 0.01 and 100, for an oxygen-rich red 
giant and a carbon star (Fig. 1). Corresponding mass-loss rates are M ~ 5x 10 -8 
to 4 x 10 -5 M yr _1 for a typical Asymptotic Giant Branch (AGB) star at solar 
metallicity, but note that the difference in luminosity between massive RSGs 
and the tip of the Red Giant Branch corresponds to a two orders of magnitude 
difference in the mass-loss rate inferred from the shape of the SED alone. 

Differences in optical depth manifest themselves most through extinction 
of stellar light at A < 3 /jm, the strength of discrete features associated with 
certain minerals mainly in the A = 10 to 30 |im region, and a general elevation 
of the spectrum over that of a naked star at A > 30 /xm. The extinction only 
becomes noticeable when ry approaches unity, as it is difficult to isolate its effect 
from the depression of the stellar continuum due to blanketing by molecular 
bands in cool giants, and from the interstellar contribution to the extinction. 
Because the absorbed photons are more energetic than the reradiated photons, 
small optical depths readily give rise to detectable amounts of excess emission at 
longer wavelengths. In particular, silicates shine efficiently around A = 10 /im, 
and the resulting emission feature can be used to detect mass-loss rates as low as 
a few 10 -8 Mq yr _1 in low-luminosity or metal-rich red giants — but here, too, 
the sensitivity is limited by molecular absorption bands and uncertainty in the 
stellar continuum (stars are not blackbodies) as well as the detailed mineralogy. 
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Figure 1. Spe ctral energy distributions compu ted with the radiation trans- 
fer code DUSTY ([Ivezic. Elitzur. fc Nenkov a 1999), for a spherically symmetric 
envelope surrounding a blackbody with a luminosity of 10 4 L Q placed at the 
distance of the Large Magellanic Cloud. Two cases are considered: T* = 3200 
K and oxygenous grains (silicates), and T* = 2800 K and carbonaceous grains 
(amorphous carbon and 5% admixture of silicon-carbide). The computations 
are performed on a logarithmic grid of optical depth from ry = 0.01 to 100. 

4. A formula for the mass-loss rate as a function of optical extinction 

For the examples shown in Fig. 1, the mass-loss rate turns out to be related to 
the visual extinction, A\ (= 1.086 Ty), as M k Ay 75 . Incidentally, replacing 
r by Ay 75 in Eq. (5) and calibrating against the computed models yields a 
practically identical constant of proportionality for both silicates and amorphous 
carbon grains. Hence: 

MfMoyr- 1 ] = 1.5 x 1(T 9 Z[Z ]-°- 5 L[L ] - 75 ^ 75 . (6) 

This provides a new, simple and versatile recipe for estimating mass-loss rates. 
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5. Relationships between mass-loss rate and infrared flux or colour 



As a shortcut to comprehensive modeling of the SED, various relationships have 
been suggested between the optical depth and a monochromatic infrared flux or 
a single infrared colour. The main reason why this is interesting, is because it 
may not be possible or economic to sample all of the SED. 

The InfraRed Astronomical Satellite (IRAS) made it possible to measure 
the long- wavelength t ail of t he du st emission from galactic red giants. This 
was hrst quantified by lJural ( 19871 ). who proposed a formula for the mass-loss 
rate as a function of the flux density at 60 fim — it depends also on the wind 
speed, distance, luminosity, dust:gas ratio, and the "mean wavelength of the 
light emerging from the star and its circumstellar dust shell" (which requires 
measurement of a si gnificant portion of the SEP). The fo rmul a was adapted to 
be use d at 25 fim bv lWhitelock. Feast, fc Catchpold (|l99ll ) and lWhitelock et al.l 
(|1994T ). Both formulae compare favourably with measurements of the mass-loss 
rate from the CO(J=1^0) rotational transition at A = 3 mm. 

Relationships between mass-loss rate and colours are appealing as they 
are distance invariant, but in the light of the previous discussion this seems a 
rather deceptive advantage as they must depend on the star's luminosity. Such 
(and other unaccounted for) dependencies are not fully appreciated when they 
are derived for relatively uniform samples of stars, for instance carbon stars 
or nearby low-mass red giants in a narrow range of metallicity. Formulae are 
considered successful if they yield mass-loss rates to an intrinsic accuracy of a 
factor two or so, but the same accuracy cannot be guaranteed if applied outside 
its validity domain (not limited by just colour range). 

I nonetheless summarise in Table 1 several formulae for the mass-loss rate 
as a function of infrared colour. Where a second reference is given the data 
presented by the former were first parameterised by the latter. Linear, square- 
root or asymptotic functions have been suggested. They share many of the 
basic assumptions, and notably they are all evaluated for stars of a (near) solar 
metallicity — even the mass-loss rates for the presumably metal-poor magellanic 
stars were derived adopting values for the dust:gas ratio and wind speed typical 
for red giants encount ered in the solar neighbourhoo d. 

The appr oaches oflLe Sidaner k Le Bertrd (|l996l ) , iLe Bertrd (ll997T).lLe Bertre fc Winters! 
(ll99ST). and oflSchutte fc Tielend (11989T) (parameterised bvl Ziilstra et alJ ll99d) 
and Groenewegen et al.l (|2007l ) (parameterised in this review) are based on mod- 
eling of the SEDs to derive mass-loss rates and then correlating with the near- 

I 1 i i 

infrared colours of the same stars. T he approach of iWhitelock et al 
(parameterised by IZijlstra et al.l Il996l ) is similar as they use the Jur 



1994) 



1983) 



formalism to derive the mass-loss rate, which is a crude way of modeling the 
SED. They found that the correlation between the mass-loss rate and 25:60 
urn flux ratio is rather poor. The mass-loss rates from the infrared data in 
Le Sidaner fc Le Bertrd (I1996T ). ILe Bertrd (119971 ). ILe Bertre Winterd (| 19981 1 
and IWhitelock et al.1 (|1994D " were all found to be in good agreement with the 
mass- loss rates from the CO(J=l — >0) line emission for the same stars. The ap- 



proach of Guandalini et all (|2006l ) (parameterised in this review) differs in that 
the mass-loss rates that were correlated with the infrared colours were not in 
themselves derived from infrared data but from the CO line emission. 



Table 1. Formulae for the 
strictions. Dependencies on 



mass-loss rate 
lum inosi ty L [L^ 



as a function of infrared colour, with a synopsis of the more pertinent re- 
are implic it unless noted othe r wise. Re ferences: l=Groenewegen et all 



(12007ft. 2=lGuandalini etafl (120061). 3^Le Bertrd (I1997T). 4=lLe Bertre fc Winters! ffl998h. 5=lLe~Sidane 



6= ISchroder. Wachter. k Winters! |2003h . 7= ISchutte fc Tielend (|1989h . 8= IWhitelock et all (|1994l ). 9= JZiilstra et al 



fc Le Bertre 




10=this review. 



jj= \ogM{x) [Mq yr 1 ] = a(x + b) c + d validity domain chemistry tp w cxp remarks ref 
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Schroder. Wachter. Wintersl (2003) compute hydrodynamical models in 



order to generate synthetic mass-loss rates and infrared colours as the model 
stars evolve. The models are computed for initial stellar masses limited to 0.63- 
1.2 Mq, and are only valid for the superwind regime of the highest mass- loss 
rates. They also predict H-K and L-M colours but these are less sensitive 
measures of the mas s-loss rate. Essentially the same models were computed by 
Winters et al.l ( 2000l ). who showed that although the evolution of mas s-loss rate 



with i nfrared colour is qualitatively comparable to that observed by iLe Bertrel 
( 19971), the results are of fset and display a different slope. 



Groenewegen ( 20061 ) computes the SEDs for a grid of carbon- and oxygen- 



rich AGB stars covering a range in mass-loss rate (similar to what was done for 
Fig. 1), and convolves these SEDs with a suite of filter curves to generate the 
expected optical and infrared photometry. If applying the appropriate scaling 
relations as discussed before, these magnitudes can be compared with actual 
measurements to estimate the mass- loss rate of a star. 



6. Complications with deriving mass-loss rates from the dust 

The description of the dust-driven wind in §2 is a highly idealized one, but often 
the observational constrains are lacking that would justify a more detailed one. 
It is thus important to be aware of the caveats, of which I list a few here. 

The dust:gas ratio is difficult to measure except for nearby objects from 
which CO line emission can be detected (realising that CO is also a minor 
constituent). In well-developed dust-driven winds, the dust condensation seems 
to be complete in the sense that all condensable material is locked up in the 
grains. This yields a dust:gas mass ratio ip ~ 0.005 for stars of an overall solar 
metallicity, with the expected scaling of Eq. (4). The validity of this notion is 
confirmed by agreement between the mass-loss rates predicted from a formula 
that was derived for a magellanic sample of oxygen-rich AGB stars and RSGs, 
and the mass- loss rates measured from the far-infr ared brightness, for a sample 
of dust-enshrouded galactic AGB stars and RSGs ( van Loon et al.l [20051 ). 



These authors also show (in their Fig. 11) that the paradigm breaks down 
in the case of giants with warm photospheres or that pulsate less vigorously. In 
the case of Betelgeuse this is known to be associated with a low dust:gas ratio 
suggesting that the dust condensation process has not reached its full potential. 
Although this regime holds vital clues to the dust formation process, it makes 
it an extremely hazardous affair to predict the dust:gas ratio and hence the 
total mass in the wind. This might gain in relative importance at low overall 
metallicity, when the carbomoxygen ratio approaches unity, or when Hot Bottom 
Burning in massive oxygen-rich AGB stars depletes both carbon and oxygen by 
burning it into the non-refractive element nitrogen. The exact composition of 
the dust and the grain size distribution then also become more uncertain and 
with it, the optical constants that determine the opacity. 

In dust-poor winds one can expect a large drift velo city between grains and 
gas, a nd Eq. (3) would need to be modified to reflect this ( Habing. Tignon. &: Tieleni 



1994). The SED shortward of its peak also depends on the acceleration of the 



wind, as well as the radial profile of the dust condensation process. This can be 
addressed in detailed non-LTE models for the line profiles of rotational transi- 
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tions at radio wavelengths ( Decin et al. 20061 ). The dusty code offers the option 
to use a hydrodynamic model for the wind, but not a radial stratification of the 
dust properties or dust:gas ratio. Stellar radial pulsations drive shocks through 
the inner parts of the envelope which cause density enhancements and variations 
in the outflow speed, further modifying the SED. Where the mass-loss rate is 
derived from the far-infrared emission from the cold dust, this can be affected 
by temporal variations in the wind speed and mass-loss rate on timescales of 
centuries to as much as 10 5 yr — comparable to the duration of the superwind 
phase or the interval between thermal pulses that occur in the nuclear burning 
layers of AGB stars. Modeling the full S EP can, in principle, recover some of 
the mass-loss history ( Groenewegen 1995), but complications with the interpre- 
tation arise if the wind has swept up a signi ficant amount of interstellar matter 
rtVillaver. Garcia-Segura. fc Manchadcl 12004 ) ■ 

The dust envelope is usually assumed to be spherically symmetric. This 
may be a good approximation for most AGB stars, but RSGs are often found to 
display deviations from this. Worse still, RSGs are encountered near their birth 
sites, and their envelopes may be emerged in a n anisotropic external irradiation 
field ( Schuster. Humphreys. &: Marengol 20061 ). Binarity is sometimes invoked 
to ex plain axi-symmetric outflows from post-AGB objects (e.g., the Red Rect- 
angle: I Waters et al.|[i~99 8). but less is known about the effects whilst still on the 
AGB. Although the geometry of axi-sym metric envelopes can be determined 
from the SED (e.g., IWhitnev et al.l 12003). the problem becomes degenerate if 
the geometry is more complex but no spa tial information is ava ilable — which 
is the norm in extragalactic objects (e.g., Ivan Loon et al.lll999T ). This may be 
especially important for low-mass red giants, which have a large peculiar veloc- 
ity with respect to the ambient interstellar medium and are th us expected to 
develop bowshocks ( Villaver. Garcia-Segura. fc Manchadoll2004l ). 



7. Mass-loss rates from atomic and molecular infrared transitions 



Molecular hydrogen has transitions at infrared wavelengths, notably at A = 28.22 
/xm (J=2— >0) with an excitation tempe rature of T Pxn = 510 K and at A = 
17.03 ftm (J=3-*l) with T exc = 1015 K (|Burton. Hollenbach. fc Tielenslll992h . 



Unfortunately, the decay times of these quadrupole transitions are in excess of 
a century and hence the expected line emission from circumstellar envelopes is 
extremely feeble. Higher-level transitions around A = 2 iim have b een used to 



study the pulsating atmospheres of cool giants ( Hinkle et al. 2000l ). but until 



the pulsation and its connection with the mass loss are better understood it is 
not possible to reliably measure the mass-loss rate from these lines. 

Absorption in the Q-branch of the acetylene (C2H2) band at A = 13.7 
/ im has been used to estimate the mass-loss rates from magellanic carbon stars 
( Matsuura et al. 20061) . in agreement with the mass- loss rates derived from mod- 



eling of the SEDs ( van Loon et al. 2006). In its current form it is simply a 



measurement of the equivalent width and hence it suffers from dependencies on 
luminosity, velocity (structure) and chemistry similar to the optical depth of the 
dust envelope. But more sophisticated modeling of the absorption band shapes 
and line shifts in high-resolution spectra of the acetylene bands at A = 3-14 /jm 
holds great promise of improvement of the method. 
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Interesting in particular for dust-poor environments, fine-structure lines of 
atomic oxygen at A = 63 jim and atomic and singly-ionized iron, silicon and 
sulphur at A = 24-36 um. can be used to estimate the gas mass in the line - 
forming region (|Haas. Glassgold. & Tielen3ll995l : lAoki. Tsuii. Qhnakal[l998h . 
But, as with the high-excitation H2 lines, a theory for the mass loss resulting 
from such an atmosphere is required in order to infer the mass-loss rate. 

An electronic transition, the triplet of atomic helium at A = 1.083 /im can be 
used to trace the bulk flow in the chromosphere that exists around warm giants 
(T* > 4600 K), much more effectively than the more comm only used optical 
lines of Ha and Can H+K (jSmith. Dupree. k, Straderj 120041 ) . Unfortunately, 
the analysis depends sensitively on the ill-defined structure of the chromosphere. 



8. Concluding remarks 

Several methods are available to determine the mass-loss rate of cool giants from 
infrared data. These can be divided in methods based on the dust extinction and 
emission, and methods based on the absorption and emission by molecules and 
atoms. The "dust" methods are ultimately linked to the shape of the spectral 
energy distribution, be it through modeling or the use of derived relationships 
between the mass-loss rate and far-infrared flux density, infrared colour or optical 
extinction. The accuracies obtained are almost never better than a factor two 
even in nearby objects with exquisite data, and can easily be worse than an order 
of magnitude in distant objects. More importantly, systematic uncertainties of 
similar magnitudes arise from lack of theory, in particular with regard to the 
structure of the molecular atmosphere and the dust condensation process. 

Acknowledgments. I would like to warmly thank the organisers and all 
participants for a very interesting and pleasant conference in the heart of Europe. 
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